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ABSTRACT 

Within the distance of 1 pc from the Galactic center (GC), more than 100 young massive 
stars have been found. The massive stars at 0.1-1 pc from the GC are located in one or two 
disks, while those within 0.1 pc from the GC, S-stars, have an isotropic distribution. How these 
stars are formed is not well understood, especially for S-stars. Here we propose that a young star 
cluster with an intermediate-mass black hole (IMBH) can form both the disks and S-stars. We 
performed a fully self-consistent ./V-body simulation of a star cluster near the GC. Stars escaped 
from the tidally disrupted star cluster were carried to the GC due to a 1:1 mean motion resonance 
with the IMBH formed in the cluster. In the final phase of the evolution, the eccentricity of the 
IMBH becomes very high. In this phase, stars carried by the 1:1 resonance with the IMBH 
were dropped from the resonance and their orbits are randomized by a chaotic Kozai mechanism. 
The mass function of these carried stars is extremely top-heavy within 10". The surface density 
distribution of young massive stars has a slope of —1.5 within 10" from the GC. The distribution 
of stars in the most central region is isotropic. These characteristics agree well with those of stars 
observed within 10" from the GC. 
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Introduction 



More than 
been found 

near-infrared observations 



100 young massive stars have 
in the Galactic centerr (GC) by 



(Krabbc et al. 199 



Paumard et~aH 120061: iLu et al l 120061: |Po et al 



2009; iBartko et al.l 120101) . These stars appear to 
reside in one or two disk s at more than 1" from 
the GC (ILu et al.1 120061: 



Paumard etHI 120061) , 



while B-type stars within 1" (S-s tars) have an 
isotro pic and thermal distr ibution ( Schodel et al.l 
iGillessen et al. 2009). Two major scenar- 
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ios have been proposed for the formation of these 
stars. One is the in-situ formation in an accretion 
disk (jLevin fc Beloborodovi 120031 ) and the other 
is the migration of a star cluster formed severa l 
parsec or more away from the GC ( GerharcfeoOll ). 

In the 1990s, only OB and Wolf-Rayet (WR) 
stars could have been observed in the GC. How- 
ever, recent advance in observation techniques 
has made it possible to identify many late-type 
stars, and now we can obtain their spacial dis- 
tribution. While the OB and WR stars have a 
power-law distribution with a slope of about —1.5, 
the slope fo r all stars including late-type giants 



is nearly too et al.l liool iBuchholz etaHl2009t 



Bartko et al.l 120101 ) . The or bits of S-stars appea r 



to be isotropic and thermal (jSchodel et al. 2003), 



but the most recent observation showed that it is 
more eccentric with th e distribution of n(e) ~ e 2 ' 6 
(Gillesscn et al . 2009). How these distributions 
have formed is the key to the understanding of the 
formation process of the young stars near the GC. 
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The mass function of the young stars have also 
been deter mined from the observations. It is quite 
top-heavv too et all EM iBuchholz etaH l2009t 



Bartko et al.ll2010t) . Formation scenarios should 



also explain the mass function. 

In-situ formation of S-stars is difficult because 
of the strong tidal field of the central supermas- 
sive black hole (SMBH). Therefore, it is neces- 
sary to carry young stars to the current location 
from somewhere outside. Recently, several scenar- 
ios have been proposed to carry stars which were 
born on a gaseous disk at around 0.1-0.5 pc from 
the GC and to randomize their orbits. The scenar- 
ios include migratio n via the gravitational to rques 
in the stellar disk (|Levinl 120071: iGrivl |2009h. mi- 



grati on as a star cluster core (jBerukoff fc Hansen 
2006), randomization of the orbital elements of S- 
aalandris fc Merritd l2009t 



Merritt et al 


2009) 


(Perets et al. 


2009h. 



or stellar-mass black holes 
and the formation of S-stars 



due t o disruptions of bin aries (|Lockmann et al 



2008; iMadigan et all 120 09). On the other hand, 



Bonnell fc Ricd (|2008l ) suggested the direct forma- 
tion of S-stars in a gaseous disk formed from a gi- 
ant molecular cloud infalling to the GC. However, 
none of them is well established. 

We have shown that star clusters can carry 
young stars to the GC, i f an I MBH forms in the 
clusters. In Fuiii et all ( 20091 ). we performed a 
fully self-consistent iV-body simulation, in which 
both the internal dynamics of a star cluster and 
the interaction between the cluster and its par- 
ent galaxy are handled correctly. The star cluster 
migrating into the GC was completely disrupted 
by the tidal force and the stars escaping from the 
cluster formed a disk structure. Before the dis- 
ruption, an IMBH formed through the runaway 
collisions of stars in the cluster. We found a new 
migration mechanism of young stars, a 1:1 mean 
motion resonance with the IMBH. The IMBH car- 
ries young stars to the GC by the 1:1 resonance 
after the disruption of the cluster. In this simu- 
lation, however, the spacial resolution around the 
SMBH was limited to 0.2 pc because of the use 
of a large softening length for the SMBH. Hence, 
it was impossible to compare the distributions of 
the young stars obtained in the simulation with 
the observed one. 

In this Letter, we report the result of a new sim- 
ulation performed using our improved code which 



does not need softening for the SMBH. In this sim- 
ulation, we can follow the orbits of stars down to 
the AU scale, where S-stars reside. We found that 
the distribution of stars obtained from the sim- 
ulation agrees very well with the observations in 
the following three points: the surface density of 
young massive stars has a slope of —1.5, the young 
stars have an extremely top-heavy mass function, 
and the orbits of young stars in the inner most 
region are thermal and isotropic. A sinking star 
cluster can explain both a young stellar disk and 
S-stars at the same time. 

We describe the method of our iV-body simula- 
tion in Section 2. In Section 3, we show the results 
of simulations. Section 4 is for summary. 

2. iV-body simulation 

We adopted a King model with the non- 
dimensional central potential of Wq — 10 as the 
Galaxy model. We plac ed a central SMBH wit h 

~ 20051) . 



a mass of 3.6 x 10 6 M(T, (lEisenhauer et al 



Our galaxy model represents the central region of 
our Galaxy. The total mass of our galaxy model 
is 5.8 x 10 7 Af Q in the real scale (excluding the 
SMBH). The number of particles is 6 x 10 6 and 
the mass of a particle is 9.7 M Q . The half-mass 
radius is 22 pc and the initial core radius before 
we place the SMBH is 0.8 pc. As a model of a star 
cluster, we also adopted a King model with the 
non-dimensional central potential Wq = 6. The 
number of particles is 64k. We started a simu- 
lation with the initial position of the cluster of 
12.5 pc from the GC and with an eccentric or- 
bit. The initial orbital velocity is 0.67 times the 
circular velocity. The total mass of the cluster is 
2.1 x 10 5 Mq and the half-mass radius is 0.16 pc. 
Initial mass function (IMF) of stars in the clusters 
is a Salpeter I MF with lower and upper cutoff at 
1 and 100 Mq (ISalpeterlll955h . We assigned each 
star a mass randomly chosen from the Salpeter 
IMF, irrespective of its position. In this simu- 
lation, we adopted collisions of stars in the star 
cluster, mass-loss from very massive stars, and 
formation of an IMBH in the clus ter. The details 
are described in (|Fuiii et alJl2009h . 

Th e calculation was performed using the Bridge 
code (|Fuiii et al.l l2007). which is a tree-direct hy- 



brid scheme. It can treat dense regions embedded 
in a larger-scale system. In this simulation, the in- 
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ternal motion of star clusters and field stars within 
0.28 pc from the central SMBH are calculated us- 
i ng a direct code, the sixt h-order Hermite scheme 
( Nitadori fc Making 2008f) . and ot her interactions 
are c alculated using a treecode (| Barnes fc Hut 
19861 ). The softening length is 3.1 x 10" 3 pc for 
field stars. We did not use softening for other in- 
teractions. We assume that stars in the distance 
200 times the Schwarzshild radius of the SMBH 
(0.14 AU) merge to the SMBH. The step size for 
the treecode is 15 years. It is short enough to re- 
solve the orbital motion at the boundary between 
the tree and direct regions. Thus, we can follow 
the orbits of stars down to the AU scale, in which 
S-stars are located. We used the opening angle 
9 = 0.75 with the center-of-mass approximation 
for the treecode. The simulation was performed 
on the Cray XT4 at National Astronomical Ob- 
servatory of Japan. We used 256 cores, and the 
total CPU time for the calculation was about 100 
days. 

3. Distribution of young stars in the 
Galactic center 

The star cluster spiraled into the GC due to the 
dynamical friction. At around 1 pc from the GC, 
the cluster was completely disrupted by the tidal 
force of the parent galaxy and the remnant stars of 
the star cluster formed a thick disk. Before the dis- 
ruption, an IMBH formed in the cluster through 
the runaway collisions of stars. A massive star 
grew via repeated collisions and finally the mass 
of the massive star reached about 16000 M Q . Af- 
ter the disruption, the IMBH continued to sink 
to the GC, and a few hundred stars which were 
the members of the cluster also migrated into the 
GC due to the 1:1 mean motion resonance of the 
I MBH. This e volut ion is the same as that shown 
in iFuiii et all (|2009l ). but in the current simulation 



line of ascending nodes at ft = 99 ° east of north. 



These values follow the result of iPaumard et al 
( 2006[ ). The arrows show the proper motions of 



we can follow the orbits of th e IMBH and stars 
down to AU scales, while in Fuiii et al. (2009) 



stars could not go below 0.2 pc because of the 
large softening length. 

Figure Q] shows the projected distribution of 
stars at the end of the run (9.82Myr). We plot- 
ted only the stars which were originally the mem- 
bers of the star cluster. The orbital plane of the 
star cluster has an inclination of i = 127° with 
respect to the plane of the sky and with a half- 



the stars. Red arrows are clockwise and green ones 
are counterclockwise. Dotted curves show the or- 
bits of stars which reached the projected distance 
less than 0.05 pc from the central SMBH. The dis- 
tribution of stars in this snapshot is very similar to 
those of S-stars, but the scale is larger in the simu- 
lation. As we stated above, these stars are carried 
by the 1:1 mean motion resonance with the IMBH 
and released from the resonance when the orbit 
of the IMBH becomes highly eccentric (see figure 
(2). Thus, the radial distribution of the carried 
stars strongly depends on the r adius at which the 



IMBH orbit becomes eccent r ic. |Baumgardt et al . 
(|2006l) . iMatsubavashi et al.l (|2007l) . and Iwasawa 



et al. (preprint) showed that the orbit of IMBHs 
becomes highly eccentric after their orbital decay 
slowed down because of the depletion of stars. The 
mechanism is as follows (Iwasawa et al, preprint): 
When the orbit of the IMBH is exactly circular, 
the orbits of stars with the semimajor axis larger 
than that of the IMBH would follow the Kozai 
cycle, which is a large-amplitude coupled oscilla- 
tion of eccentricity and inclination. In this case, 
the z-component of the angular momentum, J z , is 
conserved, since the time-averaged potential of the 
SMBH and the IMBH is axisymmetric. If the or- 
bit of the IMBH is not exactly circular, for some 
stars J 2 is no longer conserved. Thus, they can 
exchange the angular momentum with the IMBH. 
Moreover, they are likely to be ejected by the 
IMBH when they are in the prograde orbits with 
the peri-center distance comparable to the semi- 
major axis of the IMBH. On average, these stars 
carry away the orbital angular momentum of the 
IMBH and cause an increase in the eccentricity of 
the IMBH. 

We investigated the surface density profile, av- 
eraged over the last 10 snapshots with a time in- 
terval of around 8000 years. Since the young stars 
distribute mainly in the orbital plane of the star 
cluster, we plotted the projected surface density 
for inclinations i = 0°, 30°, 60°, and 90°. The 
inclination of the observed clockw ise disk at 0.1-1 
pc is 127° (jPaumard et al.ll2006l) . It is similar to 
i = 60° . Red curves in figure [3] show the surface 
density of all young stars. Green and blue curves 
show those of massive stars only with m > 5M Q 
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and to > 20 M©, respectively. The error bars 
show their standard deviations. While the mas- 
sive stars concentrate on the GC, the distribution 
of all stars is rather flat in all inclinations. The 
slope of the projected surface density profile of the 
massive stars is around —1.5 within 10", although 
it is slightly different for the different inclinations 
and the number of particles is too small to fit in 
the case of small inclinations. This result agrees 
well with the observed surface density of young 
mass i ve stars ( Bartko et al. 2010t Buchholz et al 



2009: iDo et al 



2009) 



We also investigated the mass functions of 
young stars in three intervals, r < 10", 10" < 
r < 25", and 25" < r < 35", averaged within 
the last 10 snapshots. Figure [4] shows the re- 
sult for i = 60°, which is a similar angle to the 
observed young stellar disk. The mass function 
was extremely top-heavy in the inner-most re- 
gion, r < 10", and approaches a power-law mass 
function in the outer region. These results are co n- 
sistent with the observation ( Bartko et al. 2010l) . 



As shown in iFujii et al.l (|2009h . the star cluster 
can carry massive stars efficiently to the GC be- 
cause of the following mechanisms. Massive stars 
sink to the cluster center because of the mass seg- 
regation. When a star cluster is disrupted by the 
tidal force of the parent galaxy, stars in the outer 
region of the cluster become unbound first. There- 
fore, massive stars tend to be carried close to the 
GC. After the disruption of the star cluster, some 
of the massive stars are carried further into the 
inner region by the 1:1 mean motion resonance of 
the IMBH. 

Figure O shows the probability distribution 
function of young stars within 0.3 pc from the GC 
(left panel of figure [5]) and the distribution of incli- 
nations (right panel of figured]). The distribution 
of the ecce ntricity show s a thermal distribution, 
n(e) = 2e (jHeggielll975h . T his result agrees wit h 



the simulation performed bv lMerritt et al. (2009) 



The inclination is randomized near the GC, but 
no star has an inclination larger than 90° in this 
snapshot. However, the inclination of one reso- 
nant star exceeded 90° . Several stars escaping the 
cluster by slingshot also had retrograde orbits and 
have the pericenter distance less than 0.1 pc. 



4. Summary and Discussion 

We performed a fully self-consistent A-body 
simulation of a star cluster near the GC. The star 
cluster migrated to the GC owing to the dynamical 
friction and was disrupted by the tidal force. The 
stars which are initially the members of the star 
cluster formed a disk structure. Before the dis- 
ruption, an IMBH formed in the cluster via run- 
away collisions of stars. After the disruption of 
the cluster, the IMBH continued to sink to the GC 
and stars which were caught in the 1:1 mean mo- 
tion resonance of the IMBH also sank to the GC. 
Near the GC, the spiral-in of the IMBH slowed 
because of the depletion of field stars, and the 
orbit of the IMBH became highly eccentric. At 
this stage, the stars were kicked out from the res- 
onance and the orbits were efficiently randomized 
by the non-axisymmetric perturbing potential of 
the IMBH. 

We investigated the distributions of the stars 
carried to the GC by the star cluster and the 
IMBH. We found that they agree well with the ob- 
served ones. The surface density within 10" had 
a slope of —1.5. The mass function of the young 
stars was extremely top-heavy in the inner-most 
region. The eccentricities and inclinations of the 
young stars carried near the central SMBH by the 
resonance were a thermal and isotropic distribu- 
tion, while young stars in the outer region were 
distributed in a disk. These distributions agree 
with that of S-stars and a young stellar disk. Thus, 
the distributions of "debris stars" of the sinking 
star cluster agree well with the observations. The 
star cluster scenario with an IMBH can explain 
the origin of both a young stellar disk and S-stars. 

Here, we discuss possible conditions for the for- 
mation of S-stars. In our simulation, the dis- 
tances of young stars carried by the resonance 
were larger than those of S-stars. However, how 
deeply the IMBH can carry stars to the GC de- 
pends on the distance where the orbit of the IMBH 
becomes eccentric. It occurs when the mass of 
the IMBH is c omparable to the enclosed mass of 
the field stars ( Matsubavashi et alJ l2007) and this 
distance strongly depends on the density distri- 
bution of the GC and the mass of the IMBH. 
First, we discuss the density distribution. The 
enclosed mass of the field stars in our model is 
only ~ 10% of that estimated by observed visible 
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stars (jSchodel et al.l 120071) at 0.15 pc, where the 
orbit of the IMBH became eccentric in our sim- 
ulation. If the density profile is a broken power- 
law fitted by observed visible stars , the enclosed 
mass is ~ lOOOMtn even at 0.01 pc (jSchodel et al 



20071: iGenzel et al J 120031 ). We also estimate the 
case of a power-l aw with —7/4, which is theoret- 
ically expected (jBahcall fc Wolj Il976l ). In this 
case, the enclosed mass is 8000M Q at 0.01 pc 
(2000 AU) and 1000 M at 2 x 10~ 3 pc (400 
AU). These values are smaller than the upper 
limit of the enclosed mass fro m the observatio n 
of S2, 3-4 xlO 5 M at 0.01 pc dGhez et al]l200"3) . 



Recent obse rvations show a flat density profile 
of old stars (|Po et al J 1 20091: iBuchholz et aLT f2009; 



Bartko et al 



2010). However, it is unlikely that 



this distribution reflects the true mass distribu- 
tion. There are certainly dark masses composed 
of stellar mass black holes, neutron stars, and 
white and brown dwarfs. On the other hand, 
the timescale of collisions between main-sequence 
stars is pretty small, around 0.1-1 Gyr for inner 
1 pc, if a s tellar cusp developed thro ugh thermal 
relaxation ( Duncan &: Shapiro! Il983h . Therefore, 
the lack of old stars is probably the result of col- 
lisional disruptions. Compact objects are not dis- 
rupted by collisions, and brown dwarfs are also less 
likely to be disrupted by collisions due to their 
high density. Thus, it seems natural to assume 
that the innermost region of the GC is dominated 
by dark mass. 

Next, we discuss the possible range of the 
IMBH mass. In our simulation, the IMBH is 
more massi ve than the observational upper limit, 
- 1O 4 M (jReid fc Brunthalerl [iool . but it is 
possible to form smaller IMBHs from different 
i nitial conditions or with higher mass loss rates 
(|Fuiii et al.ll2009h . liven if we assume an extreme 
mass loss rate, it is possible to form a star with a 
few thousand solar masses through runaway colli- 
sions because the star cluster near the GC is very 
compact. For the lower limit, 150 0MfT) is suffi - 



cient for the randomization of stars (jMerritt et al 
2009). Less massive IMBHs can carry stars closer 



to the GC, but it takes a longer time to mi- 
grate due to the dynamical friction. Star clus- 
ters migrate to around lpc and are disrupted 
there. Its typical timescale is 2-10 Myr for a 
star clu ster at 5-10 pc f rom the GC with around 
10 5 M,r, (|Fujii et all l2009j ). After the disruption, 



IMBHs migrate due to the dynamical friction. We 
estimated the timescale for IMBHs migrating from 
1 pc to 10~ 3 pc using the equation derived from 



Matsubavashi et alj (120071 ). We assumed that the 
Bahcall-Wolf cusp and obtained 1 Myr for the 
IMBH with 16000M©. This result is consistent 
with our simulation. We also tried broken power 
laws (jGenzel et al.ll2003l ) and obtained similar re- 
sults. The timescale is inversely proportional to 
the IMBH mass. Thus, IMBHs with a few thou- 
sand solar masses are capable for this scenario. 
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supported by Research Fellowships of JSPS for 
Young Scientist. Numerical computations were 
carried out on Cray XT4 at the Center for Com- 
putational Astrophysics, CfCA, of National As- 
tronomical Observatory of Japan. 




Fig. 1. — Projected distribution of stars at T = 
9.82 Myr. The orbital plane of the star cluster has 
an inclination of i — 127° with respect to the plane 
of the sky, with a half-line of ascending nodes at 
ft = 99° east of nor th. These value s are u sed to 
mimic the result of iPaumard et al.1 (|2006l ) . Ar- 



rows show the proper motion of stars. Red and 
white crosses show the positions of the IMBH and 
SMBH, respectively. Dotted lines show the or- 
bit of stars which reached inner 0.05 pc from the 
SMBH (projected distance). 



5 




2 4 6 8 10 



T (Myr) 

Fig. 2.— Orbital evolutions of the IMBH (red) 
and a resonant star (black). 
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Fig. 3. — Surface density of young stars as a function of projected distance from the SMBH averaged within 
the last ten snapshots. Top left, top right, bottom left, and bottom right panels show the projections with 
inclinations i = 0°,30°,60°, and 90°, respectively. Red curves includes all stars. Green and Blue curves 
include stars with m > 5M© and m > 20M©, respectively. Error bars show their standard deviations. 
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Fig. 4. — Mass functions of young stars in three 
radial intervals, r < 10", 10" < r < 25", and 
25" < r < 35", averaged within the last ten snap- 
shots. We rotated the stellar disk with i = 60° 
from the orbital plane of the IMBH for the com - 
parison with the observation ( Bartko et al. 2010f) . 
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Fig. 5. — Distribution of the orbits of young stars around the GC. Left: cumulative distribution for the 
eccentricities within 0.3 pc for the last three snapshots. Dashed curve shows a thermal distribution, n(e) = 2e. 
Right: distribution of the inclinations of young stars in the last three snapshots. 
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